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ABSTRACT 

A521 is an interacting galaxy cluster located at z — 0.247, hosting a low-frequency radio halo connected to 
an eastern radio relic. Previous Chandra observations hinted at the presence of an X-ray brightness edge at the 
position of the relic, which may be a shock front. We analyze a deep observation of A521 recently performed with 
XMM-Newton in order to probe the cluster structure up to the outermost regions covered by the radio emission. 
The cluster atmosphere exhibits various brightness and temperature anisotropies. In particular, two cluster cores 
appear to be separated by two cold fronts. We find two shock fronts, one that was suggested by Chandra and that is 
propagating to the east, and another to the southwestern cluster outskirt. The two main interacting clusters appear 
to be separated by a shock-heated region, which exhibits a spatial correlation with the radio halo. The outer edge 
of the radio relic coincides spatially with a shock front, suggesting that this shock is responsible for the generation 
of cosmic-ray electrons in the relic. The propagation direction and Mach number of the shock front derived from 
the gas density jump, M = 2.4 =t 0.2, are consistent with expectations from the radio spectral index, under the 
assumption of Fermi I acceleration mechanism. 
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1. INTRODUCTION 

Collisions between massive galaxy clusters are the most en- 
ergetic events in the present universe. Part of the kinetic energy 
released during these collisions is dissipated through supersonic 
shock fronts propagating in the intracluster medium (ICM) and 
turbulent motions. While heating the thermal component of 
the ICM, shocks and turbulence may also accelerate (or re- 
accelerate) relativistic particles (e.g., Cassano & Brunetti 2005; 
Hoeft & Briiggen 2007; Ryu et al. 2003; Pfrommer et al. 2006; 
Brunetti & Lazarian 2007; Vazza et al. 2009). Radio observa- 
tions probe these complex mechanisms through the detection of 
diffuse synchrotron emission from the ICM, in the form of giant 
radio halos, Mpc-scale radio emission in the cluster central re- 
gions, and radio relics, sharp-edged radio sources in the cluster 
periphery (e.g., Ferrari et al. 2008; Cassano 2009; Venturi 2011; 
Brunetti 2011, for recent reviews). 

Observable as sharp X-ray brightness and temperature dis- 
continuities, few shock fronts have been detected so far because 
they can only be visible in the brightest cluster regions and in 
particularly favorable projections (Markevitch et al. 2002, 2005; 
Russell et al. 2010; Finoguenov et al. 2010; Macario et al. 201 1). 
Peripheral radio relics are believed to be shock fronts that propa- 
gated far outside the X-ray bright region while still accelerating 
(or re-accelerating) electrons, which produce radio emission and 
quickly cool after the shock passes, resulting in a characteris- 
tic narrow feature (e.g., Ensslin et al. 1998; van Weeren et al. 
2011). The physics of giant radio halos is probably more com- 
plex. Radio halos plausibly result from the (re-)acceleration and 
transport of relativistic particles in large turbulent regions of 
the ICM, although many aspects of the mechanisms generating 
radio-emitting electrons remain unclear (e.g., Brunetti 2011, for 


a recent review). Sharp radio edges (and radio relics) are fre- 
quently observed at the border of giant radio halos, suggesting 
a possible link between merger shocks and the generation of 
turbulence in the ICM (e.g., Markevitch 2010; Macario et al. 
2011 ). 

A521 is a moderately distant (z = 0.247) and X-ray luminous 
(L x = (5.2+ 1.2) x 10 37 W; Amaud et al. 2000) 7 galaxy cluster, 
presenting several signatures of dynamical activity. As revealed 
already in ROSAT images, its X-ray and optical components 
appear spatially segregated, with an N-S bimodality of the 
X-ray emission (Arnaud et al. 2000) and a more complex galaxy 
number density distribution revealing two NW/SE and NE/SW 
major elongations (Ferrari et al. 2003). As further shown from 
Chandra data analysis, the ICM in A521 exhibits an irregular 
thermal structure with indications for gas heating at the interface 
between the two main gas components (Ferrari et al. 2006). 
A521 exhibits a giant radio halo that is the prototype of the class 
of ultrasteep spectrum radio halos (Brunetti et al. 2008). Since 
the halo is spatially connected to a radio relic, A521 provides 
us with an ideal test case to investigate the effects of shocks 
on the properties of the thermal and non-thermal components 
of the ICM and their connection with giant radio halos. In 
this merging cluster, a shock is suggested by the presence of an 
X-ray brightness edge on the SE side of the cluster coinciding 
with the edge of the radio relic (Giacintucci et al. 2008). The 
larger-scale cluster radio halo shows a very steep synchrotron 
spectrum supporting a picture where relativistic electrons are 
stochastically reaccelerated by the nonlinear interaction with 
turbulence in the ICM (Brunetti et al. 2008). 

7 1 0 37 W = 10 44 erg s _1 ; X-ray luminosity in the 0. 1-2.4 keV band has been 
corrected for luminosity distance assuming Hq = 70 km s -1 Mpc -1 and 
A = 0.7. 


1 


The Astrophysical Journal, 764:82 (llpp), 2013 February 10 


Bourdin et al. 


Table 1 

Effective Exposure Time of Each XMM-Newton - EPIC Observation 


XMM-Newton 

ObsIDs 

Centre Coordinates 

MOS1 Effective 
Exposure Time (ks) 

MOS2 Effective 
Exposure Time (ks) 

PN Effective 
Exposure Time (ks) 

0603890101 (S) 
0603890101 (U) 

04 h 54 m 22*00— 10° 16'30." 
04 h 54 m 22?00— 10° 16'30." 

15.7 (50.8%) 
64.4 (77.3%) 

15.7 (68.9%) 
64.5 (74.6%) 

12.8(36.1%) 
57.6 (52.5%) 


Note. The fraction of the useful exposure time after solar- flare filtering is shown in brackets. 



Arc Minutes 

Center: R.A. 04 54 19.64 Dec -10 18 02.2 
Figure 1. EPI C-XMM-Newton exposure of A521. 


The present article will focus on the analysis of a deep 
observation of A521 recently performed with XMM-Newton , 
with the particular goal of probing the ICM structure up to the 
outermost regions covered by the cluster radio halo and radio 
relic. After discussing data preparation and analyses issues in 
Sections 2 and 3, we present the various X-ray brightness and 
temperature features revealed by this observation in Section 4. 
We comment on the interplay between thermal and non-thermal 
components of the ICM in Section 5. Unless otherwise noted, 
any energy distribution is normalized as a probability density 
function, while confidence ranges on individual parameter 
estimates are 68%. In the following, intracluster distances are 
computed as angular diameter distances, assuming a A CDM 
cosmology with Hq = 70 km s -1 Mpc -1 , Qm = 0.3, and 
£2 a = 0.7. Given these assumptions, an angular separation of 
1 arcmin corresponds to a projected intracluster distance of 
232.5 kpc. 

2. OBSERVATIONS AND DATA PREPARATION 

The EPIC -XMM-Newton data set is a dual observation of 
A521, performed with focal aim point and central EPIC-MOS 
CCDs located in the SE cluster outskirts (see Figure 1). In order 
to remove the contribution of soft proton flares, we filtered the 
histogram of the photon arrival times through a temporal wavelet 
analysis. A summary of the “good” exposure time remaining on 
each of the three EPIC cameras is provided in Table 1. The 
average good exposure time is about 75 ks. 


In order to perform imaging and spatially resolved spec- 
troscopy, we binned photons in sky coordinates (k, l) and energy 
(< e ), matching the angular and spectral resolution of each focal 
instrument. To map the surface brightness of extended sources, 
these photon counts may have to be normalized for spatial and 
spectral variations of the telescope effective area and detector 
exposure times. We thus associate an “effective exposure” ar- 
ray, E(k, l, e ), with the photon event cube. Expressible, e.g., in 
s deg 2 , E(k, l, e) is computed as a linear combination of CCD 
exposure times, tcco(k, /, /?), related to individual observations 
p, with local corrections for useful CCD areas, accD^k, /, p), 
Reflexion Grating Spectrometer (RGS) transmissions, 8 tr^Qs, 
and mirror vignetting factors accvik, /, p) 9 

K 

E(k , l,e) = ^2 fccD (k, /, p) x Aa mmor (k, /, e, p) 
p = i 

X t v rgs (k , /, e, p) 
x a C cD(k , /, p). ( 1 ) 

3. DATA ANALYSIS 
3.1. Background Noise Modeling 

The cluster emissivity must be separated from an additive, 
spatially extended and mostly stationary background noise, 
including false photon detections due to charged particle- 
induced and out-of-time events, but also the cosmic X-ray 
background (CXB), and some Galactic foreground components. 

The XMM- EPIC background is dominated by the particle 
component, which is modeled from observations performed in 
the filter wheel closed (FWC) position during revolutions 230 
to 2027 for the EPIC-MOS cameras, and 355 to 1905 for the 
EPIC-PN camera. Following an approach proposed in, e.g., 
Kuntz & Snowden (2008) or Leccardi & Molendi (2008), this 
model sums a quiescent continuum to a set of florescence emis- 
sion lines convolved with the energy response of each detector. 
It is completed with a residual emission associated with soft 
protons, presently only detectable in the case of the EPIC-PN 
camera and modeled as a power-low spectrum normalized to 
1.4 cts deg -2 s -1 in the 0.5-1. keV band. To account for the two 
different spectral shapes in the soft and hard bands, the quies- 
cent continuum is modeled as the product of a power law with 
an inverted error function that increases in the soft band. We set 
the emission line energies to the values reported in Leccardi & 
Molendi (2008), while the soft proton residual is modeled using 
an additional power law. Presumably due to differences in the 
collecting areas of the imaging and readout detector regions, 
the EPIC-MOS quiescent continuum exhibits a small emissiv- 
ity gradient along the RAWY CCD coordinate, which has been 


8 EPIC-MOS detectors share a common optical path with the RGS. 

9 Information about these instrumental effects has been obtained from the 
XMM-Newton-EPIC Current Calibration Files (CCFs). 
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measured and taken into account in the model. Because the flu- 
orescence lines exhibit a more complex spatial variation (Lumb 
et al. 2002; Kuntz & Snowden 2008), we modeled the emissivity 
distribution of the most prominent lines 10 from the wavelet fil- 
tering of a set of FWC event images in narrow energy intervals 
around each line. 

Secondary background components include the CXB and 
Galactic foregrounds. Since they are associated with real pho- 
ton detections, these components are corrected for the effective 
exposure. The CXB is modeled with an absorbed power law of 
index y = 1.42 (see, e.g., Lumb et al. 2002), while the Galactic 
foregrounds are modeled by the sum of two absorbed thermal 
components accounting for the Galactic transabsorption emis- 
sion (TAE; kTl = 0.099 keV and kT2 = 0.248 keV, see Kuntz 
& Snowden 2000). We estimate the emissivities of each of these 
components from a “joint-fit” of all of the background compo- 
nents in a region of the field of view located beyond the boundary 
of X-ray emission in the SE cluster outskirt, but covered by the 
central MOS CCDs (see also Figure 2). This estimates yields 
13.4, 28.0, and 29.5 cts m -2 deg -2 s _1 in the 0.5-1. keV band 
for the two transabsorption and CXB components, respectively 
(X 2 /dof= 1.23). Our background model, Ab C k(&, /) Fb C k 
(k, l, e ), eventually includes a contribution for the EPIC-PN 
out-of-time count rate, which is estimated in each energy band 
as 6.3% of all photon counts registered along the CCD columns. 

3.2. Spectroscopic and Surface Brightness Measurements 

To estimate average ICM temperatures, kT, and metal abun- 
dances, Z, along the line of sight and for a given location of 
the field of view (k, /), we add a source emission spectrum to 
the background model and fit the spectral shape of the result- 
ing function, A ev t(&, l)F QWt (kT, Z, Ah, e ), to the photon energy 
distribution registered in the energy band (0.3-12 keV): 

A evt (k, l)F QYt (kT, Z, A h , e) 

= E(k, /, e) x NicM(k, l)FicM(kT, Z, Ah, e) 

+ A/bck^bck(fc> l, e). (2) 

In this modeling, the source emission spectrum F\ C M(kT, Z, 
Ah, e) assumes a redshifted and Ah absorbed emission modeled 
from the Astrophysical Plasma Emission Code (APEC; Smith 
et al. 2001), with the element abundances of Grevesse & Sauval 
(1998) and the neutral hydrogen absorption cross sections of 
Balucinska-Church & McCammon (1992). The Ah value has 
been fixed to 4.9 x 10 24 m -2 from measurements obtained 
near A521 in the Leiden/Argentine/Bonn Survey of Galactic 
Hi (Kalberla et al. 2005). It is altered by the mirror effective 
areas, filter transmissions, and detector quantum efficiency (see 
footnote 8), and convolved by a local energy response matrix 
M(k,l, e, e') computed from response matrixes files (RMF) 
tabulated in detector coordinates in the XMM-Newton-EPIC 
calibration data base. * 11 

To compute images and radial profiles of the intracluster 
gas distribution, we estimate a cluster surface brightness map, 
X x (k, /), from photon counts registered in a given energy band 
and corrected for effective exposure and additive background. 
Assuming an average ICM energy distribution, (Ccm(^)) ? we 


10 Namely, the Al, Si and Cu, Ni complexes as for the EPIC-MOS and 
EPIC-PN cameras, respectively. 

1 1 EPIC response matrixes are computed from canned RMFs corresponding to 
the observation period provided by the XMM-Newton Science Operation 
Centre. 



1 10 
Energy (keV) 



Figure 2. Background spectrum observable in the A521 outskirts. Light blue: 
particle background. Cyan blue: cosmic X-ray background emission. Blue and 
violet: TAE emission (kTl = 0.099 keV and kT2 = 0.248 keV, see Kuntz & 
Snowden 2000, and details in Section 3.1). Green: residual soft proton emission. 
Red and black: overall fit and data set. 

(A color version of this figure is available in the online journal.) 

define X x (k, /) as a function of an effective exposure map, 
e(k,l) = Z e {Ficu(e))E(k,l,ey. 

_ , . Ee »evt (k, l, e) - ./VbckCE l ) Xe E bc k(E /, g) 

xC ’ ’ ~ e(k. 1) ' ( 

All parameters of (Fi C m(<0) = Fi C M(kT 0 , Z 0 , V H , 0 , e) are 
practically determined from spectral fitting of the main cluster 
emission spectrum: kT 0 = 6.7 keV, Z 0 = .4Z 0 , Ah, 0 = 
4.9 x 10 24 m -2 , while X x (k, l) is estimated in a “soft” energy 
band ([.5-2.5] keV 7 ) in order to lower the dependence of 
(Eicm(^)) on kT 0 . 

3.3. Surface Brightness, ICM Density, 
and Temperature Profiles 

In the following, surface brightness and temperature profiles 
have been extracted within cluster sectors oriented approxi- 
mately along the surface brightness gradients. We derived the 
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radial surface brightness X x (r) by averaging the surface bright- 
ness X x (k, /) of Equation (3) in each profile annulus composed 
of N pixels ( k , /), as follows: 


£ x (r) : 


k,l 


(4) 


Since the background contribution N^(k, 0 ^bck (k, l, e) 
is estimated within a much larger area of the field of view 
than any sector annulus used to derive X x (r), we neglected any 
systematic uncertainty related to its modeling and estimated the 
variance on Z x (r) from a weighted mean of the local Poisson 
fluctuations in X x (k, /): 


, x2 _ 1 d^evt ’ 0 2 _ 1 ^evt(£> 0 

CTZ, r ~~ Nj^~e(kjW~ - 


(5) 


Projected temperatures £T(r) and the associated confidence 
interval 8kT(r) have been computed within each annulus by 
fitting a uniform emission model to the data set. To do so, 
we averaged the emission models of Equation (2) associated 
with each pixel (k, l) of the annulus, and estimated the model 
parameters kT(r ), Z(r), and A^r), via a x 2 minimization. 

These brightness and temperature profiles have been used 
to model the underlying density and temperature of the ICM, 
assuming spherical symmetry of the cluster atmosphere in 
the vicinity of the features of interest. This was undertaken 
by projecting and fitting parametric distributions of the three- 
dimensional (3D) emission measure 12 , n p n e , and temperature, 
T(r), to the observed profiles. In this modeling, projected bright- 
ness profiles are convolved with the XMM-Newton focus point- 
spread function (PSF), while projected temperatures are com- 
puted assuming the “spectroscopic-like” weighting scheme pro- 
posed in Mazzotta et al. (2004). In Section 4.3, the ICM emission 
measure (see footnote 12) and temperature profiles across two 
shock fronts have been modeled by step-like distributions with 
a common jump radius ry 


{ D Z n 2 o( r / r j) 2m ’ r < r j 
\n\(r t r jY 2 " 2 , r > r,- 


T{r) — 


[D^To, r < r j 

I T„, r > rj 


3.4. Imaging and Spectral Imaging 
3.4.1. Imaging 


( 6 ) 


(7) 


An image of the cluster is presented in the top panels of 
Figure 3. To obtain this image, X x (k, /), we corrected the EPIC- 
XMM-Newton raw photon image for spatially variable effective 
area and background flux, following Equation (3). The point- 
like sources have also been modeled by means of an isotropic 
undecimated B3-spline wavelet analysis (see, e.g., Starck et al. 
2007) and subtracted from the image. 

A map of anisotropic details in the ICM structure is shown 
in the bottom-left panel of Figure 3. To create this image, we 
subtracted a wavelet filtered map of the photon rate, E x (k, /), 
from the photon rate itself, then smoothed the residual image 
with a Gaussian function of typical width FWHM = 20 arcsec. 


12 More precisely, the ICM emission measure per volume unit. 


The wavelet filtering was performed by means of a soft 3 a 
thresholding of B3-spline wavelet coefficients, with the signifi- 
cance thresholds being directly computed from the raw — Pois- 
son distributed — photon map, following the multiscale variance 
stabilization scheme introduced in Zhang et al. (2008). 

3.4.2. Spectral Imaging 

In order to map the ICM temperature in A5 21, we used the 
EPI C-XMM-Newton data set and applied the spectral-imaging 
algorithm detailed in Bourdin et al. (2004) and Bourdin & 
Mazzotta (2008, hereafter B08). Following this algorithm, a 
set of temperature arrays kT(k, /, a) with associated fluctuations 
< 7kt(k , l, a) are first computed on various analysis scales a , then 
convolved by complementary high-pass and low-pass analysis 
filters in order to derive wavelet coefficients. The wavelet 
coefficients are subsequently thresholded according to a given 
confidence level in order to restore a denoised temperature map. 
Here, the signal analyses have been performed over six dyadic 
scales within an angular resolution range of 8a = [1.7-110] 
arcsec. This was undertaken by averaging the emission modeled 
by Equation (2) within overlapping meta-pixels (k,l,a) and 
computing the kT(k, l, a) and cr^(k, /, a) arrays by means of a 
likelihood maximization. The resulting ICM temperature map 
shown in Figure 3 was then obtained from a B 2-spline wavelet 
analysis (see B08 for details) with coefficients thresholded to 
the lcr confidence level. 

4. ICM THERMODYNAMICS 
4.1. Intracluster Gas Brightness and Thermal Structure 

The X-ray photon image of Figure 3 reveals the complex 
morphology of the intracluster gas in A521. On large scales, a 
northern sub-cluster with a comet shape is apparently falling 
onto the main component. The photon image also reveals 
the strongly irregular morphology of the surface brightness, 
presenting various edges indicated with dashed lines. Some 
of these brightness jumps have been enhanced in the bottom- 
left image of anisotropic details. They are also noticeable on 
the surface brightness profiles of Figures 4 and 5. At the 
interface between the two main interacting cluster components, 
we observe two bow-shaped brightness jumps, CFi and CF 2 , 
joining with each other to form a low brightness cross-shaped 
feature. A third brightness jump with a higher curvature radius, 
Si, is crossing the southern cluster component from SE to NW, 
while a fourth one, S 2 , is visible at the southeast cluster outskirts. 

The ICM temperature map of Figure 3 is strongly irregular, 
and presents various noticeable features. The northern sub- 
cluster is clearly cool (kT ~ 4.5 keV). The interacting region 
separating this cool core from the main cluster to the south 
appears hotter (kT > 7 keV) and strongly disturbed. The cross- 
shaped brightness depression observable on the photon image 
seems to coincide with a hot cross (kT ~ 9 keV), in particular, 
along the brightness jump CFi. The southern part of the main 
cluster is cooler (kT ~ 4 keV) than the interacting region, in 
particular, to the south of the brightness jump Si . 

Bringing together the brightness and temperature maps of 
Figure 3, we observe that the two brightness jumps CFi and 
CF 2 are associated with temperature increases as the brightness 
decreases, while the brightness jumps Si and S 2 are associated 
with a temperature decrease. The two jumps CFi and CF 2 are 
thus likely to be cold fronts separating the densest parts of the 
two sub-clusters from their interacting region, while Si and 
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Figure 3. EPI C-XMM-Newton observation of A521. Top panels: photon rate image in the 0.5-2. 5 keV band. Photon counts in this image have been corrected for 
spatially variable effective area, background flux, and wavelet-detected point-like sources. Bottom-left panel: anisotropic details in the ICM emissivity map. These 
details have been enhanced from subtraction of a wavelet denoised map to the photon rate (further details are provided in Section 3.4.1) Bottom-right panel: ICM 
temperature map obtained from wavelet spectral imaging. Prominent brightness jumps are indicated by dashed lines on the photon rate image and temperature map. 
(A color version of this figure is available in the online journal.) 
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Figure 4. Projected gas brightness measured across two cluster sectors intercepting the brightness jumps CFi and CF 2 . The projection of a step-like gas density 
distribution (Equation (6)) convolved with the XMM-Newton PSF is superimposed as a dotted line, assuming density jump amplitudes of 1.7 ± 0.1 and 2.1 ±0.1 for 
CFi and CF 2 , respectively. 
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Figure 5. Top panels: projected gas brightness and temperature profiles measured across two cluster sectors intercepting the brightness jumps Si and S 2 , as shown in 
Figure 7. Bottom panels: ICM density and temperature profiles modeled as step-like 3D distributions matching the projected profiles (see also Equations (6) and (7)). 
Dispersions on these profiles have been estimated from random realizations of the data set and corresponding models, each profile envelope delimiting 68% of the 
realizations with closest x 2 distance from the original data set. The projection of these distributions is reported as a dotted line on the projected profiles. 



S 2 are probably shock fronts propagating outward from the 
colliding clusters. 

4.2. Cold Fronts and Shock Heating at the Interface Between 
Two Interacting Sub-clusters 

One of the most striking features seen in our data is a 
cross-shaped brightness depletion separating the two colliding 
sub-clusters. This feature also corresponds to a temperature 
and entropy enhancement, in particular just outside the two 
cold fronts, CFi and CF 2 . What we see is probably shocked 
gas with high entropy being squeezed by the converging cool 
core remnants and flowing around the densest part of the two 
interacting clusters without penetrating the two cold fronts. The 
projected layer of shocked gas would thus exhibit maximal 
temperature and entropy near the two cold fronts, where it 
is tangentially intercepted by the line of sight. This shocked 
gas layer might also partly overlay in projection the main sub- 
cluster from its boundary delineated by cold front CF 2 to the 
southern brightness jump, Si. A possible interpretation for the 
origin of this hot gas flow is illustrated in Figure 6. Originally 
located at the cluster boundary (if there is one), the high entropy 
gas may have been shock heated between the two clusters that 
were starting to interact. It would then expand over the cluster 
atmosphere, following the shock fronts presently propagating 
outside the cluster cores. One of these shock fronts might be 
observed to the south of the main cluster as the brightness and 
temperature jump, Si. 

4.3. Shock Fronts Propagation to the Cluster Outskirts 

The 2D gas brightness and temperature maps of Figure 3 
suggest that the brightness jumps Si and S 2 are shock fronts 
propagating outward from the cluster center. Located at various 
distances from the cluster center, these two shock fronts might 
have been developed during two successive cluster collisions. 
In order to analyze ICM thermodynamics across these jumps, 
we extracted the brightness and temperature profiles shown in 
Figure 5, corresponding to the two sectors of Figure 7. An 
estimate of the cluster photon counts in each sector region is 
provided in Table 2. The brightness jumps Si and S 2 exhibit the 
typical shape of a projected spherical density jump convolved 
with the XMM-Newton PSF. We model the underlying gas 
density and temperature profiles as two step-like functions with 
a common jump location, following Equations (6) and (7). A 


Cluster core 

Cluster boundary 

Cold front 

Shock front 

Shock front (projection) 

Shocked gas 

Figure 6. Tentative interpretation of the ICM thermal and entropy structure 
observed in the central region of A521 . Left: early stage of a two-cluster merger: 
the cluster boundaries start to collide and develop two shock fronts propagating 
within the densest regions of each cluster. Meanwhile, the clusters develop two 
cold fronts while pushing the higher entropy gas away from their interacting 
region. Right: the shock fronts have now propagated to the most external regions 
of the interacting clusters, but could not penetrate the two cool cores. A shocked 
gas region with high entropy remains at the interface between the two cold 
fronts. 

(A color version of this figure is available in the online journal.) 


discussion about the validity of the assumption of the ICM 
spherical symmetry in the vicinity of the shock fronts is provided 
in the Appendix. The gas density and temperature distributions 
corresponding to this model are reported under the projected 
profiles in Figure 5. The 3D density and temperature jumps 
associated with these distributions are reported in Table 3, with 
confidence intervals estimated at 68% of a parameter sample 
matching several random realizations of the data set. 

The direction of the temperature jumps is consistent with the 
shock front interpretation. The cold front hypothesis would in- 
stead imply a temperature increase across the jumps (D T < 1), 
which is excluded by the data. Assuming two shocks propagat- 
ing outward in the main cluster, one should be able to estimate 
the shock Mach numbers from the Hugoniot-Rankine density, 
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Table 2 

Estimated Cluster Photon Counts Within the 0.3-5. 0 keV Band in the Regions Shown in Figure 5 


Detector 


Sector 1 



Sector 2 



Region 1 

Region 2 

Region 3 

Region 4 

Region 1 

Region 2 

Region 3 

EPIC-MOS 1 
EPIC-MOS2 
EPIC-PN 

985(79.4%) 

1016(80.4%) 

2086(74.0%) 

518(86.1%) 

552(87.6%) 

1083(82.3%) 

592(88.9%) 

578(89.4%) 

1146(86.4%) 

879(90.4%) 

868(91.0%) 

1557(87.3%) 

140(19.6%) 

87(14.9%) 

192(13.6%) 

173(37.0%) 

144(36.4%) 

284(29.2%) 

232(49.8%) 

232(49.8%) 

420(42.9%) 


Note. The fraction of the total counts is shown in brackets. 


Table 3 

Density, Temperature Jumps, and Mach Numbers Estimated Across the Shock Fronts Si and S 2 


Shock Front 

Density Estimates 


Temperature Estimates 


Jump Amplitude Mach Number 

Jump Amplitude 

Mach Number 

51 

5 2 

1 48+0.02 
A * 4 °-0.11 
9 64+0-13 
Z - D4 -0.15 

1 to+0.02 

1 — 0.08 

2 42 +019 
Z * 4Z — 0.19 

1 70 + 0 . 63 
A ' ' — 0.38 
4 47+B.O6 
1 -3.00 

1 76 +0 - 49 
A ' 7 —0.35 

3 - 40 -l 6 9 9 2 


1000 


500 


o 

£ 0 


-500 


-1000 


Figure 7. Photon rate image of A521 extracted in the 0.5-2. 5 keV band. The 
image has been re-binned to a 6.8 arcsec angular resolution in order to enhance 
the brightness jumps Si and S 2 . The two annular sectors show the two cluster 
regions where the temperature and brightness profiles of Figure 5 have been 
extracted. 

(A color version of this figure is available in the online journal.) 



-4 -2 0 2 4 

Arc minutes 

R. A. (0) 04 54 1 4.41 Dec(0) -1 0 1 5 02.4 


temperature, or pressure jump conditions across the fronts. Such 
Mach number values are reported in Table 3. The Mach num- 
bers independently estimated from the density and temperature 
jumps are consistent with each other, though estimates from the 
temperature jumps have larger uncertainties. We will hereafter 
use Mach number estimates for both shocks S 1 and S 2 , from their 
density jumps: M s \ fP = 1. 33^08 anc * M$ 2 ,p = 2.42 zb 0.19. 


5. NON-THERMAL ICM EMISSION 

A521 hosts a radio relic in its southeastern peripheral region 
and a rare low-frequency giant radio halo. In order to investigate 
the interplay between thermal and non-thermal components 


of the ICM emission, the 240 MHz radio image obtained 
from observations performed at the Giant Metrewave Radio 
Telescope has been superimposed on the X-ray photon and ICM 
temperature maps of Figure 8. 

5.1. The A521 Radio Halo 

The A521 radio halo was discovered from low-frequency 
observations at the GMRT (240, 325, 610 MHz; Brunetti et al. 
2008, see also Figure 8) and then studied in detail through a deep 
follow-up Very Large Array observation at 1 .4 GHz (Dallacasa 
et al. 2009). Its very steep spectrum, with spectral index a ~ 1 .9 
between 325 and 1400 MHz, suggests magnetohydrodynamic 
turbulence to be responsible for the in situ re- acceleration of 
the relativistic electrons (Brunetti et al. 2008). The radio halo 
covers the cluster central region, exhibiting an EW elongation 
and reaching the radio relic to the southeast. When excluding 
the relic region, the halo appears to be spatially correlated with 
the cluster X-ray emission. There is an even better correlation 
between the radio brightness and the hottest regions of the 
ICM — in particular, the radio brightness exhibits a quick drop 
across the Si shock. 

The complex thermodynamics of the ICM in the cluster 
center hint at the possible origin of the turbulence that may 
re-accelerate non-thermal particles in the halo. The two cold 
fronts CFi and CF 2 may have developed K-H instabilities at 
large angles from the main cluster collision axis. As suggested 
by the spatial correlation between shock-heated regions and 
the radio emission, turbulence may alternatively have been 
generated behind the two shocks Si and S 2 , now propagating 
to the cluster outskirts. In addition, the merger disturbance 
has likely generated turbulence within the two sub-cluster core 
remnants. 

5.2. Shock Front Propagation and the Radio Relic 

A521 has been known to host an SE radio relic observed 
at various frequencies (Ferrari et al. 2006; Giacintucci et al. 
2006, 2008, hereafter, G08). As shown in G08, the integrated 
synchrotron radiation in the relic exhibits a power-law spec- 
trum with spectral index a — 1.5 in the frequency range 
235-4890 Mhz, with evidence of a steepening of the radio spec- 
trum with increasing distance from the eastern edge. As further 
noted in GO 8, the outer edge of the radio relic coincides with the 
X-ray edge S 2 , which we have shown in this work to be a shock 
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Figure 8. Left panel: anisotropic details in the X-ray emissivity (same as Figure 3) overlaid with emissivity iso-contours in the 240 MHz radio band (Giant Metrewave 
Radio Telescope; Brunetti et al. 2008). Right panel: ICM temperature map overlaid with the same radio emissivity iso-contours as for the left figure. 

(A color version of this figure is available in the online journal.) 


front propagating towards the cluster outskirts. As observed in 
several peripherical radio sources of galaxy clusters (see, e.g., 
Briiggen et al. 2012 for a recent review), these facts support 
the shock electron (re-)acceleration being at least partly respon- 
sible for the radio emission from the relic. Assuming diffuse 
shock acceleration for the origin of the emitting electrons, in 
the test particle approach, the slope of the injection spectrum of 
cosmic rays is related to the shock Mach number, M, by (B land- 
ford & Eichler 1987) 8 t = 2 (M 2 + 1 )/(M 2 - 1). This leads to 
a spectrum of electrons in the downstream region with slope 
8 = 8i + 1 (implying a synchrotron spectral index a = 8 t / 2), 
taking into account radiative losses and assuming stationary 
conditions. G08 thus predicted a shock propagation with Mach 
number M ~ 2.3 from their measurement of a. The steepening 
of the radio spectrum with increasing distance from the eastern 
edge further allowed them to predict a shock propagation to the 
cluster outskirts. The propagation direction and Mach number 
of the shock front S 2 , M p = 2.4 zb 0.2 (see Section 4.3), are 
fully consistent with this hypothesis. 

Furthermore, note that the X-ray edge corresponding to the 
shock front seems to extend in the north-south direction more 
than the radio relic. A first interpretation for this limited extent 
of the relic might be that the shock would re-accelerate pre- 
existing relativistic electrons in the ICM. In this case, the radio 
relic could reflect the spatial and energy distribution of the 
pre-existing electrons across the shock front. In line with this 
hypothesis, recent analyses (Kang & Jones 2007; Kang & Ryu 
20 1 1 a, 20 1 1 b) suggest that the presence of pre-existing particles , 
in addition to the thermal pool, can significantly increase the 
average efficiency of the particle acceleration and the expected 
synchrotron emission at weak shocks (M ^ 3). Differences 
in extension between the shock and the radio relic might 
alternatively indicate some changes in the efficiency of electron 
acceleration changes along the shock front, possibly due to local 
variations of the Mach number (e.g., Hoeft et al. 2008). In 
this respect, the radio relic in A521 is located at the extremity 



Arc minutes 

R. A. (0) 04 54 11.40 Dec(0) -10 14 47.0 

Figure 9. Projected galaxy density distribution derived from photometric ob- 
servations performed at the CFH telescope (dark color indicate higher densities, 
Dressier algorithm; see Ferrari et al. 2003). Pink iso-contours: radio emis- 
sivity in the 610 MHz band (Giacintucci et al. 2008). Black iso-contours: 
X-ray emissivity in the 0.5-2.5 keV band (curvelet denoising of the 
XMM -Newton image). 

(A color version of this figure is available in the online journal.) 


of the NW/SE major galaxy alignment evidenced in Ferrari 
et al. (2003, see also Figure 9), where, indeed, recent accretion 
of sub-cluster material may have produced inhomogeneities in 
the ICM. 
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A521 is a complex cluster system where optical analyses have 
revealed at least three galaxy groups to the SE and four groups 
to the NW, including the cluster BCG group coinciding with 
the X-ray peak (Ferrari et al. 2003). The X-ray morphology of 
the BCG group suggests an infall along an NNW-SSE direction 
(projected onto the sky plane), which is slightly offset with 
respects to the major NW-SE galaxy alignment (Ferrari et al. 
2006). The cluster atmosphere exhibits various brightness and 
temperature edges associated with cold fronts and shock fronts, 
which were revealed by our XMM-Newton data. 

The two main interacting gas components in the central 
region of this system are separated by a region of gas with 
lower density, higher temperature, and entropy. We interpret 
this feature as a flow of high-entropy gas being squeezed by two 
converging sub-cluster cores that are delimited by cold fronts. 
We suggest that this high entropy gas was heated by shocks 
formerly developed when the two gas components started to 
interact. One of these shocks is currently observed to the south 
of the main component, with Mach number M p = 1.33+ qq 2 8 . 
The hot gas region separating the two interacting components 
appears spatially correlated with the cluster radio halo. The 
development of turbulence in the hot gas flowing between the 
two cool cores may be responsible for high-energy electron re- 
acceleration, yielding the radio halo emission. Merger shock 
propagation and/or cold fronts may have contributed to the 
development of these instabilities. 

A shock front is observed at the southeast cluster outskirt. 
An X-ray brightness edge there was hinted at by Chandra data 
(G08), though the statistical significance was marginal. The 
orientation of this shock front and its large distance from the 
cluster center suggest that it is associated with a cluster collision 
that has occurred prior to the current two-component interaction. 
Our Mach number for this shock, M p = 2.4 zb 0.2, is consistent 
with that expected from the spectrum of the radio relic in G08 
under the assumption of the Fermi I acceleration mechanism. 
As observed in X-ray follow-ups of other radio relics — A3667 
(Finoguenov et al. 2010) and RXCJ1314.4— 2515 (Mazzotta 
et al. 2011) — its detection supports the shock electron (re)- 
acceleration to be at least partly responsible for the radio 
emission from the relic. The detection of a polarization of 
the relic would be an additional support for this process, 
complementary to the extension of its synchrotron spectrum 
to very high radio frequencies and to evidence for spectral 
steepening downstream of the shock (GO 8). Delimited by the 
shock front, the radio relic seems, however, to subtend only 
a fraction of the shock front. Differences in the spatial extent 
of a radio relic and its companion shock front have also been 
observed in the colliding cluster RXCJ1314.4— 2515 (Mazzotta 
et al. 201 1), where a radio relic seems to be confined to a small 
section of the shock front presumably distorted by a nonuniform 
gas flow. These differences may thus reflect variations of the 
efficiency of particle acceleration across the shock that could 
be driven by local variations of the Mach number and shock 
velocity. Deeper X-ray or SZ observations may enable us to 
investigate this hypothesis, though the present XMM-Newton 
image does not seem to provide evidence of any strong variation 
in the amplitude of the surface brightness edge, and thus in the 
shock Mach number. An alternative hypothesis is that the radio 
relic would reveal local inhomogeneities in the properties of the 
pre-existing relativistic electrons that would be re- accelerated by 
the shock passage. The observed connection between the radio 


halo and the relic may suggests that pre-existing relativistic 
electrons have first been accelerated by turbulent gas motions 
responsible for the radio halo emission, then re- accelerated at 
the shock front. 
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APPENDIX 

APPENDIX MATERIAL 

The ICM density and temperature distributions intercepting 
the shocks Si and S 2 have been modeled in Section 3.3 as two 
step-like functions, assuming the shock center of curvature and 
the ICM centroid to coincide with each other. 

The X-ray image of Figure 3 seems, however, to show that the 
shock fronts Si and S2 are less curved than the closest cluster 
brightness isophotes. To investigate the systematic uncertainties 
inherent to our spherical symmetry approximation, we alterna- 
tively tried to model the shock front and the ICM density as two 
spherical distributions with distinct centers. Assuming these two 
centers to be located in the plane of the sky, the ICM emission 
measure is now expressed per volume unit, as 


\n p n e ](r) 


D Wo (£) ” ’ Y < ''J 

n o (£) " ’ r > r J 


(Al) 


where r and r' refer to the norm of each radius vector in the 
shock and ICM frames, respectively. Introducing d and d ' , 
the projection of these radius vectors onto the sky plane, a 
surface brightness profile intercepting the shock is obtained by 
the integration of Equation (Al) along the line of sight: 


Z(jc) = 2 x 


/*0max pCXD PX+8x 

7<£ min J 0 Jx 


A (T(d, l))[n p n 


tfpmin JO Jx—Sx 

x (J, d', /, <p)dd dl d(j>. 


p'^e J 


(A2) 


where d f is related to J as a function of r Q , the distance separating 
the shock from the center of the ICM distribution and 0, 
the angle separating the projected radius vector to the shock 
propagation axis (« d ' = y/{rj — r 0 ) 2 — 2 (rj — r Q )d cos(0) + J 2 , 
see also Figure 10). In addition to an ICM density slope, v, and 
the shock curvature radius, and density and temperature jumps, 
rj, D n , and D T , respectively, the ICM emission measure thus 
depends on an asphericity parameter: 1 — r 0 /rj. 

We tried to invert [n p n e ](r) and its parameters from a min- 
imization of the x 2 distance separating X(x) (Equation (A2)) 
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Figure 10. ICM geometry across the shocks. Left: cluster volume cut along the line of sight. Right: cluster volume cut within the sky plane. 


Table 4 

ICM Asphericity Parameters Across the Shock Fronts S 1 and S2 


Shock 

ICM Asphericity 
(1 — r Q /r j) 

ICM Slope 

G) 

Shock Curvature 

(r,0 

Density Jump 

(At) 

Temperature Jump 
(D T ) 

Mach Number 
(Derived from D n ) 

Si 

0 oo +0 0 ° 
u - uu -o.oo 

0 50 +0 0 ° 
u,:>u -o.oo 

0 45+0-06 
0.18 

0 70 +0,29 

u - 7 -0.50 

0 51+ 016 

U,:)i -0.34 

0 1 4+0- 12 
U>i ^-0.02 

1 C9+O.OI 
A * JZ -0.01 
1 ^9+0.01 
A * JZ -0.01 
1 ^9+0.01 
A * JZ -0.01 

1 40+0.14 
1 — 0.07 

1 46+0.15 
A — 0.07 

1 C4+0.01 

^^-o.oe 

1 sn+0-70 
A *° U -0.38 

1 9 ^r+ 0.62 

A " ' — 0.36 

1 90+O.66 

t./O— o.36 

1 26 +0 ’ 10 
1,z — 0.05 

1 24 +0 - u 
A ‘ z ^— 0.05 

1 49+O.OO 
1 -0.05 

S 2 

0 oo +00 ° 
U>UU -0.00 

0 50 +0 0 ° 
u ' :,u -o.oo 

0 5 1 +0- 02 
' J - DL - o.oi 

0 oo +0 - 26 
U,UU -0.00 

0 oo +01 ° 
U,UU -0.00 

0 00 +013 
u - uu -o.oo 

9 9C+0.02 
z - z - 0.00 
9 9C+O.O2 
z ' z - 0.01 

9 26+0- 01 
Z ' Z — 0.01 

991 +0.14 
Z ' ' —0.12 
9 94+0.25 
z * ' — 0.19 
9 74+ 0.11 
Z - /4 -0.14 

4 - 50 -2 5 9°2 

4-36 7 2 3 8 6 2 

4 Q5+4-45 
■3.50-241 

2 51+0.22 
Z '^ — 0.16 
2 54+0.41 
z ‘ J ^-0.25 
2 55+0.17 
z *^ — 0.20 


from the X-ray surface brightness profiles extracted across each 
shock front (see Figure 5). Since some of the searched pa- 
rameters degenerate with one another, we first performed this 
inversion by fixing the asphericity parameter to 0 and 0.5, cor- 
responding to shocks located at distances of rj and 2 x r ; from 
the cluster center, respectively. We subsequently left all param- 
eters free to vary and report the results of our measurements in 
Table 4, with the confidence interval on each parameter being 
estimated from the 68 % of a parameter sample matching sev- 
eral random realizations of the data set. As expected, the shock 
curvature radius, density, and temperature jumps obtained when 
fixing the asphericity to zero are consistent with their estimates 
derived from the spherical model of Section 3.3. A marginal 
difference in the amplitude of the density jump is still notice- 
able, since Equation (Al) yields D n = 1.3 W 0.07 f° r shock Si, 
while Equation (6) yields D n = 1.33+qq 2 8 . This difference is 
probably related to the lack of any variation of the ICM density 
slope across the shock, following Equation (Al). Fixing the as- 
phericity to 0.5 instead of 0 also marginally affects the density 
jump, essentially due to the degeneracy between the ICM as- 
phericity and density slope. This degeneracy is noticeable in the 
case of Si, the shock front observed with the highest statistics. 
Leaving the ICM asphericity free to vary yields estimates of 
0.45 and 0.75 in the cases of Si and S 2 , respectively, consistent 
with the shock curvature radii observed on the X-ray image of 
Figure 3. The shock density, temperature jump, and Mach num- 
bers derived from these various assumptions are, in any case, 
consistent with one another, and with their estimates obtained 
from the spherical model of Section 3.3. Given the limited statis- 


tics available, it is difficult to break the degeneracy between the 
ICM asphericity, ICM density slope, and shock curvature radius 
in the vicinity of the shocks. For simplicity, we consequently 
adopted the spherical model of Section 3.3 in order to derive the 
amplitudes of the density jumps and Mach numbers of the two 
shocks Si and S 2 . 
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